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ABSTRACT 

Context. X-ray observations of unabsorbed active galactic nuclei provide an opportunity to explore the innermost regions of super- 
massive black hole accretion discs. 

Aims. Our goal in this paper is to investigate the central environment of a Seyfert 1.5 galaxy IRAS 05078+1626. 
Methods. We studied the time-averaged spectrum obtained with the EPIC and RGS instruments onboard XMM-Newton. 
Results. A power law continuum (photon index T ^ 1.75) dominates the 2-10 keV energy range. A narrow iron Ka spectral line is 
clearly seen, presumably originating in a distant torus, but no broad relativistic component was detected. However, the power law and 
the iron K a line alone do not provide a satisfactory fit in the soft X-ray band whose spectrum can be explained by the combination 
of three components: a) a cold photoelectric absorber with column density N H ~ 10 21 cm~ 2 , This gas could be located either in outer 
parts of the accretion disc, at the rim of the torus or farther out in the host galaxy; b) a warm absorber with high ionization parameter 
(logf ~ 2.2) and column density (N H ~ 10 24 cirT 2 ); c) an ionised reflection where the reflecting gas could be either in the inner wall 
of a warm absorber cone or in an ionised accretion disc. 

Conclusions. The first X-ray spectroscopic measurement of IRAS05078+1626 unveils some of the standard ingredients in Seyfert 
galaxies, such as a power law primary continuum, modified by reflection from the accretion disc and by the effect of complex, multi- 
phase obscuration. However, data constrains the accretion disc, if present, not to extend closer than to 60 gravitational radii from the 
black hole. 
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IRAS 05078+1626 is a nearby Seyfert 1.5 galaxy. Before its 
identification as an infrared source it was also known as CSV 
6150 (Catalogue of Suspected Variables). Its position on the sky 
is / = 186.1 and b = -13.5 in the Galactic coor dinates. The 
cosmological redshift of this galaxy is z ~ 0.018 dTakata et alJ 
11994 . It had never been spectroscopically examined in X- 
ray prior to the observation discussed in this paper. However, 
it was detected in a number of X-ray surveys, such as the 
all-sky monitoring o f the INTEGRAL IBIS/ISGR I instrument 
dSazonov et al.ll2007l). the SWIFT BAT instrument dAiello et alJ 
I2008t iTuelleretaU 120081) . and the RXTE Slew Survey (XSS) 
dRevnivtsev et al.ll2004 

The X-ray spectroscopic properties of intermediate Seyferts 
are rather elusive: both obscured Type 1 and unobscured 
Type 2 active gala c tic nuclei (AGN) hav e been reported (e.g., 
Capp i et all 120061: iBrightman & Nandral 120081: iBianchi et ail 
2009bl) . It has been suggested that intermediate Seyfert galaxies 



are seen at intermediate inclination angles between pure "face- 
on" Seyfert Is and pure "edge-on" Seyfert 2s, which follows 
directly from the orientation-based AGN un i fication scenarios 
dAntonucci & Miller! 1 1985b lAntonuccilll993l: lUrrv & Padovanil 
1995). For this reason, X-ray spectroscopy of type 1.5 Seyferts 
may provide clues to the nature and geometrical distributions of 
optically thick gas surrounding the active nucleus, the latter be- 
ing the fundamental ingredient behind the unification scenarios. 
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IRAS 05078+1626 is included i n the FERO project 
("Finding Extreme Relativistic Objects"; Longi notti et alj [2008) 
with the aim of establishing the fraction of a relativistically 
broadened Ka iron lines in the spectrum of a complete flux- 
limited sample (2-10keV flux > 1 mCrab). 

This paper is organised as follows. In section lXTI we describe 
new observation by XMM-Newton and the corresponding data 
reduction. In sections 12.21 - 12.51 we present the results derived 
from the X-ray spectra in ~ 1-10 keV energy range. The results 
are discussed in section[3]and the conclusions are summarised in 
section|4] 



2. Results from XMM-Newton observation in 2007 

2.1. Observations and data reduction 

The XMM-Newton observation of IRAS 05078+1626 was per- 
formed between 2007 August 21 UT 22:24:49 and 22 UT 
15:3 5:43 (Obs. #050209 0501). T he EPN and both M PS cam- 
eras dStriider et al.1 d2001l) for PN; iTurner et all (1200 ll) for MOS) 
were operating in the smal l window mode. The RGS cam- 
eras dden Herder et al.ll200ll) were operating in spectroscopic 
mode. The spectra were r educed with the SAS software version 
9.0.0 (Ga briel et al.ll2.004h . Intervals of high particle background 
were removed by applying count rate thresholds on the field-of- 
view (EPIC, single events) and CCD #9 (RGS) light curves of 
0.35 cts/s for the PN, 0.5 cts/s for the MOS and 0. 15 cts/s for the 
RGS. The exposure time after data screening is w 56 ks for MOS, 
« 40 ks for PN and « 58 ks for RGS, respectively. The patterns 
0-12 were used for both MOS cameras, and patterns 0^1 (i.e. 
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single and double events) for the PN camera. The source spec- 
tra were extracted from a circle of 40 arcsec in radius defined 
around the centroid position with the background taken from an 
offset position close to the source. The two MOS spectra and the 
related response files were joined into a single spectrum and re- 
sponse matrix. Finally, the PN and MOS spectra were rebinned 
in order to have at least 25 counts per bin and to oversample 
the energy resolution of the instrument maximally by a factor of 
three, while the RGS spectra were left unbinned. Consequently, 
different statistics were used in fitting the spectra - the traditional 
y 1 statistics to fit the PN and MOS spectra and the C-statistics 
dCashl 19761) for all fits i ncluding RGS data. For the spectral anal- 
ysis, we used XSPEC (lArnaudl 19961) version 12.5, which is part 
of the HEASOFT software package version 6.601 
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2.2. Timing properties 

The PN light curve of the source is shown in Fig. Q] We have 
divided the energy range into two bands and checked the light 
curve behaviour in each of them, as well as a hardness ratio, 
which we defined as the ratio of the counts at 2-10keV to the 
counts at 0.2-2keV. The energy ranges were chosen for sam- 
pling different spectral components, as indicated by the energy 
where the continuum starts deviating from a power law model 
that describes the hard X-ray spectrum (see Sect. 12. 31 ). The hard- 
ness ratio stays almost constant during the observation, suggest- 
ing that no significant spectral variations occur, although the 
source flux increased by around 20%. 
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Fig.l. EPIC-PN light curves in the 0.2-2 keV band (upper 
panel) and 2-10keV band (middle). The hardness ratio HR is 
defined as the ratio of the counts at 2-10keV to the counts at 
0.2-2 keV and presented as a function of time in the lower panel. 
The bin time is as 2048 s. 

To confirm this conclusion also for narrow spectral features, 
such as the iron emission line, we compared the PN spectra ex- 
tracted during the first and the second halves of the observation 
(see Fig. [2j. The spectra correspond to the lower/higher source 
flux because the flux is increasing nearly monotonically during 
the observation. We calculated the ratio values of the two data 
sets and fit them with a simple function f(E) = aE + b using the 
least square method. The fitting results are a — -0.004 + 0.003 
and b = 1.12 ± 0.01 with the sum of the residuals^- 2 = 111 for 



Fig. 2. PN spectrum extracted from the first half of the observa- 
tion with the lower source flux (blue) and from the second half 
of the observation with the higher source flux (black). The ratio 
of the two spectra is presented in the lower panel. 
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Fig. 3. Total PN spectrum with the background level showing 
that the signal-to-noise ratio is very good up to high energies. 



190 degrees of freedom. When we set a = 0, the fitting results 
are comparably good with b — 1.11 ±0.01 and^ 2 = 113. The ra- 
tio of the spectra is plotted in the lower panel of Fig. [2] Because 
no significant spectral differences are evident, we analyse the 
time-averaged spectra hereafter. 

2.3. Mean spectral properties 

The signal-to-noise ratio is very good up to high energies 
(Fig. [3J, so we fit the EPIC spectra spectra in the full energy 
range where they are well calibrated (0.35-12 keV). The X-ray 
continuum is described by a power law model at energies above 
2 keV, although the iron line at E = 6.4 keV is present (Fig. |4j». 
The photon index of the power law is Y =* 1.49(1)0 In this 
and all subsequent models we included absorption by Galactic 
gas matter along the line of sight with column density «h = 
0.188xl0 22 cm 2 . T his value is from the L eiden/ Argentine/Bonn 
HI measurements (Kal berla et alj 120051) . We used the TBABS 



http://heasarc.gsfc.nasa.gov 



2 All presented errors represent 90 % confidence level for a single 
interesting parameter, and the errors quoted in brackets are related to 
the last digit in the number. 
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Fig. 5. Ratios of the simple power law model (the same as in Fig. |4j to the data in different energy bands: left: at lower energies, 
right: in the iron line band, where the narrow Ka line at the rest energy E = 6.4 keV is prominent (observed at E — 6.29 keV due 
to the cosmological redshift). Black crosses correspond to the PN data points, while the red crosses correspond to the MOS data 
points. 
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Fig. 4. XMM-Newton PN (black) and joint MOS (red) spectrum 
of IRAS05078 + 1626 described by a simple power law model 
absorbed by Galactic neutral hydrogen in the line of sight with 
«h = 0.188 x 10 22 crrT 2 . The photon index of the power law is 
T = 1.49. The model reveals an apparent excess at E — 6.4 keV 
associated with the iron line Ka and some wiggle-like residuals 
at lower energies. A more detailed view of the data residuals in 
these parts of the spectrum is shown in Fig. [5] 



sorbers in Seyfert galaxies. The spectrum can be also affected by 
the so-called soft excess, which can be either ca used by reflec- 
tion b y the ionised surface of the accretion disc (Crumm y et al.l 
2006) or by partially ion ised and Doppler smeared absorption 
dGierliriski & Donell2006l) . 

The spectral residuals reveal certain discrepancies between 
the PN and MOS spectra (see Fig. [5] with the data/model ra- 
tios of both spectra, PN and MOS, with the identical model pa- 
rameters of the spectra). The level of discrepancy is, however, 
comparable to the level of systematic u ncertainties in the cross - 
calibration between the EPIC cameras (StuhlingeretaD|2006). 
Nonetheless, we conservatively analyse the EPIC spectra sepa- 
rately. We use the same models for both spectra but allow the 
values of the model parameters to be different. The values of the 
photon index using the simple power law model differ from each 
other when fitting the spectra independently, resulting in a harder 
PN spectrum with T — 1.60(1) compared to the MOS spectrum 
with T = 1.54(1), ignoring the energies below 2keV and also 
between 5.5-7.5 keV. Although the absolute value of these spec- 
tral index measurements does not have a direct physical mean- 
ing, given the simplicity of the model applied on a small energy 
band, the comparison between them is illustrative of the quality 
of the cross-calibration between the EPIC cameras. Differences 
of the order of AF 0.06 in the hard X-ray band are cons istent 
with current systematic uncertainties (Stuhlinger et al. 200j|). 



model dWilms et al]l2000l) to fit the absorption produced by the 
Galactic interstellar matter. 

We applied the simple TBABS*POWERLAW model to both 
PN and MOS spectra. The x 2 value is 3557 with 528 degrees 
of freedom (x 2 /v = 6.7) in the 0.35 - 12.0 keV energy range. 
The spectra differ from the power law model not only around 
E = 6.4 keV but also at lower energy band 0.35 -2.0 keV (adding 
a Gaussian line model to fit the iron line improves the fit only 
to x 2 / v - 3384/524 = 6.5). Residuals against this model are 
shown in Fig. [5] 

The residuals at lower energies are usually attributed to a 
warm absorber, i.e. absorption by tota ll y or partially ionise d 



matter; see, e. 



iiy_ . 

Netzer et al.l d2003l) : iBlustin et all <|2005); 



iRrongold et al.l (120071) for more information about warm ab 



2.4. RGS spectrum 

We jointly fit the unbinned first-order spectra of the two RGS 
cameras with the same model's parameter values except the 
overall normalisations. The continuum is well-fitted by the sim- 
ple power law model with the photon index F = 1.57. We 
searched further for narrow emission and absorption lines in the 
spectrum using several ZGAUSS models with the intrinsic width 
cr set to zero. We calculated the 90 % confidence interval for a 
blind search, as P = Po/Mrfai. where A^triai = AWs/3 = 3400/3 
and P Q = 0.1. For the RGS data P = 8.8 x 10~ 5 , to which 
AC = 22.4 corresponds assuming the Student probability dis- 
tribution. The only line fulfilling this criterion by improving 
the fit about AC = 31.7 is an emission line at the energy 
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Fig. 6. Residuals of the PN data from the 'baseline' model in- 
cluding both ionised reflection and absorption (lower), only re- 
flection (middle), and only absorption (upper). 



E = 0.561 ±0.001 keV (22. 10 ± 0.04 A) and the equivalent width 
eV. We identify it with the forbidden line of the O VII triplet 
(£ LA b =0.561 keV). 

2.5. EPIC spectrum 

The forbidden line of the OVII triplet is clear signature of a pho- 
toionised plasma. No significant features were detected that may 
be expected alongside the O VII (f) line, if it were produced in 
a collisionally ionised plasma, such as the resonance line in the 
OVII triplet or the OVIII Lya. This led us to try and explain the 
residuals against a power law model in the soft X-rays as effect 
of intervening io nised absorption gas. We us ed the XSTAR model 
version 2.11n7c dKallman & Bautista 2001m1 to calculate a grid 
of tabular models with the input parameters constrained from 
the preliminary data analysis with simple models whenever pos- 
sible (photon index T « 1.7, density p < 10 14 cirT 3 , luminosity 
L < 10 44 ergs~', column density 10 19 cirT 2 < tin < 10 25 cm~ 2 , 
and ionization parameter -5 < log£ < 5). 

A single-zone warm absorber component modifying the 
power law continuum dramatically improved the fit from^/v = 
1850/270 = 6.9 to ^ 2 /v = 402/270 = 1.5 for the PN spec- 
trum. The ionization parameter converged to a very low value, 
and we found that this almost neutral absorption can be success- 
fully reproduced with ZPHABS, which is a simpler model than 
XSTAR, so we preferred this possibility. The addition of another 
warm absorber zone improves the fit to x 2 l v = 320/266 = 1.2 
for the PN spectrum, and it requires the ionization parameter 
log£ = 3.9. We checked that adding another warm absorber 
zones does not improve the fit significantly. The residuals from 
the model (see Fig. El upper panel) reveal an extra emission that 
remains at low energies, as well as around the iron K a line band. 

These features can come from reflection of the primary ra- 
diation on the surf ace of the accretion disc, so we added the 
REFLION model (iRoss & Fabianl l2005h . which calculates the 



3 http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 



Energy [keV] 

Fig. 7. The 'final' model. The total model is shown in black 
(solid line), the primary radiation is red (dashed), the REFLION 
components are blue (dotted) for cold reflection and magenta 
(dot-dashed) for ionised reflection. 



ionised reflection for an optically thick atmosphere with constant 
density. We examined the significance of the addition of the re- 
flection component into the complex model of the PN spectrum 
by the statistical F-test. The low value of the F-test probability 
(5 x 10~ 15 ) strongly favours this additional model component. 
The best fit was now^ 2 /v = 246/264 = 0.95 for the PN spec- 
trum. We hereafter call this model the 'baseline' model; in the 
XSPEC notation: tbabs x zphabs N x xstar x (powerlaw 

+ REFLION + ZGAUSS). 

The parameter values of the 'baseline' model are presented 
in the Table Q] The quoted errors of the parameters represent 
a 90 % confidence level for a single interesting parameter. The 
measurement is obviously affected by a much larger systematic 
error, which, however, could be properly quantified only if we 
knew the "right" model. The value of the power law photon in- 
dex increased toF* 1.8 compared to the simple model applied 
to the data in Sect. 12.31 because we included of the additional 
local absorption in the model. The data residuals from the model 
are shown in the lower panel of Fig. [6] In the same figure, we 
also show residuals from best fit performed with the 'baseline' 
model, excluding the ionised absorption (middle panel) and the 
ionised reflection component (upper panel). 

The narrow iron K a line with the rest energy E = 6.40 + 
0.01 keV, the width cr = 0.06+0.03 keV, and the equivalent width 
EW = 82+15 eV evidently represents cold reflection. This sug- 
gests an origin of this spectral component in the outer part of 
the disc, or from the torus. The cold reflection is also supposed 
to contribute to the soft part of the spectrum with the individual 
emission lines. For this reason, we replaced the Gaussian pro- 
file in the 'baseline' model with another REFLION component 
(called as REFLION 2 in the Table [TJ with the same values for 
the photon index and abundances as the REFLION 1 model com- 
ponent. The ionization parameter was kept free during the fit- 
ting procedure, but it very quickly converged to its lowest value 
£ = 30 (log£ = 1.477). The advantage of the REFLION model 
compared to the other available reflection models is that it also 
includes the soft X-ray lines, with the disadvantage in this case 
that the ionization parameter cannot be set to zero. 

This 'double reflection' model, in the XSPEC notation 

TBABS X ZPHABS N X XSTAR X (POWERLAW + REFLION + RE- 
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Table 1. Parameters of the 'baseline' and 'final' models. 
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Model 
component 


Model 
parameter 


'baselim 
PN 


j' model 
MOS 


'final' 

PN 


('double reflectio 
MOS 


n') model 
PN+MOS+RGS 


ZPHABS 


;i H [10 22 cm" 2 ] 


0.104+S 


f) 19Q+0.007 
-0.007 


0.102!°;°°' 


120 + °°° 8 

-0.005 


0.106!°™ 


XSTAR 


ft H [10 22 cm" 2 ] 
He/He Solar - Ca/Ca Sokr 

Fe/Fesolar - Ni/Ni S olar 


130! 2 ° 

9 o+0.1 
-0.1 

i(f) 

2 +01 

"■^■-0.2 


1 70 +2 ° 

9 4+0.1 
-0.1 

Kf) 


120! 3 ° 

9 9+1.4 

Kf) 

1 2 + ° 3 


150 +7 ° 

1 J "-20 
Kf) 

9 +0 - 2 

"■^-0.2 


130!?° 

2-5! ;° 
Kf) 
1 1 + ° 2 


POWERLAW 


r 

normalisation 


i o 1 +0.03 
1 - 81 -0.05 

(7 ± 1) x 10- 4 


(7 + 1) x 10- 4 


1 7C+0.10 

1 - ' -0.03 

(6± l)x 10- 4 


1 74+0.07 
(7 + 2) X 10- 4 


1 76+ 004 

1 - '"-0.02 


REFLION 


r 

logf 

Fe/Fesoiar - Ni/Nisoiar 
normalisation 


1.81 (b) 

3 +0 - 2 

-0.2 

0.2 (b) 
(3 ± 2) x 10-" 


1.80 (b) 

i 9+O.2 

-0.2 

0.1 (b) 

(3 ± 2) x 10-" 


1.75 (b) 

3 0+° 2 

J -"-0.2 

1.2 (b) 
(2± 1) x 10-' 


1.74 (b) 

3 +2 
0.9 (b) 
(1 + 1) x 10- 9 


1.76 (b) 

3 + ° 1 

J -"-0.2 

1.1 (b) 


REFLION 2 


normalisation 






3 ± 1 x 10- 7 


4 + 1 x 10- 7 




ZGAUSS 


E [keV] 
<x [keV] 

z 

normalisation 


6 40 + ° 01 

"■^"-0.01 

06 + ° 03 

"■"-0.04 

0.018 (f) 
(3 ± 1) x 10" 5 


6 44+ - 02 

"• H ^-0.02 

02 +0 04 

"•" -0.02 

0.018 (f) 
(3 + 1) x 10- 5 


- 


- 


- 


2 / 

x h 




246/264 


405/243 


256/266 


404/244 


C/v= 1551/1347 



Note: The sign (f) after a value means that the value was fixed during the fitting procedure. The sign (b) means that the parameter value was bound 
to the value of the corresponding parameter of the previous model component. The sign "-" means that the model component is not included in 
the total model, while dots in the right column only mean that there are more values related to the individual spectra which are not necessary to be 
all shown in the table. 



FLION), does not significantly improve the fit goodness over the 
'baseline' model (with^ 2 /v = 256/265 = 0.96 for the PN spec- 
trum), but it represents a more self-consistent astrophysical pic- 
ture. Therefore, we call the 'double reflection' model as 'final' 
model. In contrast to the 'baseline' model, it does not require 
subsolar iron abundances, see the Table [TJ where the parameter 
values for this model are presented. The 'final' model with each 
component separately drawn is shown in Fig. [7] All the plotted 
components are absorbed by a warm absorber surrounding the 
central accretion disc and two kinds of cold absorber - one from 
Galactic interstellar matter and one from local absorber in the 
host galaxy. 

Some model parameters were not allowed to vary during the 
fitting procedure. The redshift of the ionised absorber was fixed 
to the source cosmological value, because leaving it free yields 
a negligible improvement in the quality of the fit. Second, we 
used the same iron abundances across all the components in the 
model. 

In the 'final' model, the warm absorber ionization parameter 
is consistent with the ionised reflection component. This result 
is also presented in Fig. [8] where the contour lines related to the 
lcr, 2cr, and 3<x levels of^ 2 between the ionization parameters 
of the two model components are present. 

Table [2] summarises flux values of the individual compo- 
nents of the 'final' model for both PN and MOS spectra for 



Table 2. Flux values of the 'final' model and its individual com- 
ponents. 



Model 
component 


Flux at 0.5 -2 keV 

[10- 12 ergcm- 2 s-'] 
PN MOS 


Flux at 2- lOkeV 

[10-' 2 ergcnT 2 s~'] 
PN MOS 


total model 
unabsorbed model 
POWERLAW 
REFLIONj on 
REFLION cold 


7.05!° : ° 3 7.oo-;:«; 

16.6!° 2 16.7!° 2 
13 6 + ° 1 14 3 +01 

ij'"-0.2 l*t.j_ j 

1 q+0.2 1 1+0.1 
-0.2 -0.1 

1 1+0.1 1 1+0.1 
^ -0.1 1 -0.1 


25 +al 25 4 +01 
J -0.2 J ' -0.2 

25 7 + °- 2 26 5 + ° 2 

-0.2 -0.2 

22.6!°;?, 23.8!°;) 
1 6 + °' 9 + ° 1 

x -"-0.2 "- 7 -0.1 

1 c+0.1 1 O+0.2 
1.3_ 01 i.o_ 01 


^cold * 


0.12!°°: 0.07!°;°| 
0.07!°;°; 0.08!°,;°! 


06 + ° 01 03 + ° 01 
"■""-0.01 "•" -0.01 

06 + ° 01 07 + ° 01 
"■""-0.01 "-"'-o.oi 



* the ratios of the reflection component flux values to the flux value of 
the total unabsorbed model (sum of the primary and reflected radiation). 



two energy bands, 0.5 - 2keV and 2 - lOkeV, and also shows 
fractions of the reflection radiation to the total emission (sum 
of the primary and reprocessed radiation). The flux ratio is al- 
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1.5 2 2.5 3 3.5 4 4.5 

lo 9 ^XSTAR 



Fig. 8. The contour plot of the ionization parameters of the RE- 
FLION model, representing the ionised accretion disc, and of 
the XSTAR model, representing the warm absorber in the 'fi- 
nal' model. The individual curves correspond to the 1 cr, 2 cr, 
and 3 cr levels. The position of the minimal value of y 2 found by 
the fitting procedure is marked by a cross. The corresponding^ 2 
values are given in the plot. 
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Fig. 10. The best-fit values of x 1 statistics for the inner disc ra- 
dius parameter, which we obtained by gradually stepping it from 
the horizon radius to the outer radius of the disc (40QR e ). The 
dashed line is the 90% confidence level for one interesting pa- 
rameter. 
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Fig. 9. The joint fit of all spectra of the XMM-Newton instru- 
ments - PN (black), MOS (red), RGS 1 (magenta), and RGS 2 
(blue), together with the model residuals. 



most equally shared between the cold and ionised reflection 
components, and its value is in total R < 0.2. The absorption- 
corrected luminosity values of the source in the same energy 
bands are L 0.5-2 keV = (1.21±0.02)xl0 43 ergs _1 andL2-iokeV = 
(1.87 + 0.02) x 10 43 erg s _1 , respectively. 

We also used the 'final' model for a joint fitting of all the 
XMM-Newton instruments (PN, MOS and both RGS spectra) 
together. The parameter values were bound among all the spec- 
tra, only normalisation factors were allowed to vary. The good- 
ness of the joint fit is given in C-statistics because the RGS 
data are unbinned (and each individual bin contains only a few 
counts). The result is C — 1551 for a number of degrees of free- 
dom v = 1347. All the spectra, together with the residuals, are 
shown in Fig. |9j and the corresponding parameters in the last 
column of Table [1] 



3. Discussion 

3.1. Constraints on the location of the absorbers 

In this section we discuss a possible location of the ab- 
sorber's system in the 'final' model. Photoelectri c absorption 
is almost invariably observ e d in Type 2 Seyfert s dAwaki et al.l 
1 199 It iTurner et al.l Il997t iRisaliti et al.l 120021) and gener- 
ally attributed to an optically thick matter responsible for 
orie ntation-dependent classification in AGN un ification scenar- 
ios (lAntonucci & Millerill985l:lAntonuccilll993h . 

Because the IRAS 05078 + 1626 galaxy is probably viewed 
under an intermediate inclination between unobscured Seyfert 
Is and obscured Seyfert 2s, the torus rim may also intercept the 
line of sight to the AGN and absorb the radiation coming from 
the centre. The cold absorption can, how ever, also be associate d 
with the interstellar matter of the galaxy dLamastra et al.l l2006). 

Both reflection components are inside the ionised absorber 
in the 'final' model. The geometrical interpretation is that the 
cold reflection occurs on the outer parts of the disc or the inner 
wall of the torus. Reflection on the nearer part of the torus is 
heavily absorbed by the torus itself, so only radiation reflected 
on the farther peripheral part of the torus can reach the observer 
after passing through the warm absorber. However, an alterna- 
tive scenario, in which the cold reflection is unaffected by the 
warm absorber, i.e., TBABS x ZPHABS x [REFLlON co i d + XS- 
TAR x (POWERLAW + REFLlON disc )], is also acceptable with 
X 2 /v = 265/265. 

The lack of constraints on the va riability in the warm ab- 
sorbed features dKrongold et alJl2007h . caused by the moderate 
dynamical range of the primary continuum, as well as statistical 
limitations in our spectra, prevents us from precisely constrain- 
ing the location of the warm absorber. 

3.2. Constraints on the location of the ionised reflector 

The ionised reflection might occur either at the inner wall of a 
warm absorber cone or on the accretion disc. Even in the lat- 
ter case, the reflection cannot occur arbitrarily close to the black 
hole. In this section we investigate the constraints of the accre- 
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tion disc location and structure, which can be drawn from the 
lack of the significant relativistic blurring of the disc reflection 
component. 

We convolved the ionised r eflection component with the 
fully relativistic KYCONV model dDovciak et alJl2004h . Two as- 
sumptions about the disc emissivity were considered. First, the 
radial part of the intensity decreases with the power of the disc 
radius q (I cc r~ q ), where the value of q was allowed to vary be- 
tween 2 and 3.5. Second, the angular dependence was assumed 
to be isotropic which seems to be appropriate app roximation for 
our si tuation of an X-ray irradiated accretion disc ( Svobod a~et alJ 
l2009h . 

We examined the expected confidence levels of the best-fit 
values of the disc's inner radius by stepping this parameter in 
the whole range of its possible values - from the horizon to 
the outer disc radius, which we set to 400 gravitational radii 
(fl g = GM/c 2 ). The results are shown in Fig. [TU] At the 90 % 
confidence level, the accretion disc is not allowed to extend 
closer to the black hole than 60R g . 

The "relativistic blurring method" would be less appropriate 
in looking for the imprints of the innermost parts of the accre- 
tion disc if the disc were too highly ionised (l og £ ~ 4) and 
the n arrow reflection features were not present (Ross & Fabian 
120051) . However, the ionization parameter value of the reflection 
component is not so high in the 'final' model and the dominant 
feature is the intermediately ionised iron line (E as 6.7 keV). If 
we assumed a stratified disc with the ionization state decreas- 
ing with the radius from the centre, the hydrogen-like iron line 
would be also expected to appear in the spectrum (as an interme- 
diate stage between the over-ionised and mildly ionised contri- 
bution). Because it is not detected in the data, the accretion disc 
truncation provides a more reasonable explanation of missing 
signatures of the relativistic blurring. 

3.3. Mass accretion rate 

Disc truncation is expected in low-luminosit y AGN where 
the inner accretion flow is advection-dominated dNaravan & Yil 
Il994t lEsin et all 1 1991 iNaravan & McClintockl 120081 and ref- 
erences therein). The transition from the outer standard ac- 
cretion disc may occur, e . g., via the d i sc eva poration mech- 
anism dMever et al l 120001: Liu & Taaml 120091) . The observa- 
tional evidence of a truncate d accretion d isc in l ow-luminosity 
AGN was reported e.g. by lLasota et al.l (119961) : lOuataert et al.l 
(119991) . However, its prese nce is also suggested in some obser- 
vations of Seyfert galaxies dLu & Wandl2000j:lDone et alJl2000t 
I Chian g & Blae sl 120031; iM arkowitz & Reevesf l2009h and even a 
quasar dMatt et al.l20 05) where the luminosity value is estimated 
as a half of the Eddington value. Generally, it is expected that 
the lower the luminosit y, L/Lem, the larger transition radius (see 
lYuan & Naravanl2004 and references therein). Furthermore, we 
investigate whether the disc truncation hypothesis is consistent 
with the IRAS 05078+1626 luminosity. To have these quantities 
in Eddington units, we first estimated the mass of the black hole. 

IRAS 05078+1626 belongs to the sample of the infrared- 
selecte d Seyfert 1 .5 gala xies observed by a 2.16 m optical tele- 
scope dWang et alj |2006) where, among others, the velocity dis- 
persion in the O III emission line was measured. The correlation 
between the O III line width and the mass of active galact ic nu- 
cleus was discussed in iNelsonl (120001) and IBorosonl ([2003). The 
value from the optical measurements, cr 0III a 130km s -1 , cor- 
re sponds to the m ass M ~ 4 x 1O 7 M using a correlation plot 
in IBorosonl d2003l) . The scatter of the correlation is somewhat 
large with the reported limit of a factor of 5 for an uncertainty in 



the black hole mass determination, so the value only provides an 
order of magnitude estimation. 

The value of the Eddington luminosity is Lndd - 1.3 x 
10 38 M/M o ergs _1 w 5 x 10 45 ergs _1 for the given value 
of the mass. We used luminosity-dependent corrections by 
iMarconi et al.l (120041) to estimate the bolometric luminosity of 
IRAS 05078+1626 from the X-ray luminosity. Its value is L « 
5 x 10 44 ergs~' « lO^'Ledd- Correspondingly, the mass accre- 
tion rate, M - L/c 2 , is sub-Eddington with M » 0.1M E dd- This 
val ue is typical of less lu minous Seyfert galaxies (see for exam- 
ple [Bianchretay[2009 a), and is consistent with the disc trunca- 
tion hypothesis. 

4. Conclusions 

The X-ray continuum spectrum of the Seyfert galaxy IRAS 
05078+1626 is dominated by a power law with a standard value 
of the photon index (F = 1.75 in the 'final model'). The resid- 
uals from the power law continuum can be interpreted in terms 
of the warm absorber surrounding the accretion disc, and the re- 
flection of the primary radiation from the ionised matter and on 
the cold torus. The outgoing radiation is absorbed by cold matter 
(«h ~ 1 X 10 21 cm 2 ), which can be either located in the inner 
side of the torus or caused by gas in the host galaxy. The type 
of the galaxy determined from the previous infrared and optical 
research is Seyfert 1.5, suggesting that the active nucleus could 
be seen at large inclination, consistent with either interpretation 
or even allowing a combination of both. 

The ionised warm absorber occurs in the central part of the 
AGN. Its column density was found to be «h t 1 x 10 24 cm~ 2 , 
which is a rather high value compare d to the warm abs orbers 
detected in the other Seyfert galaxies (Blu stin et al.l l2005). This 
may be because we are looking through a longer optical path of 
a conical non-relativistic outflow due to the high inclination of 
the system. The ionization parameter of the warm absorber is 
log^wA = 2.5 + 1.0, which is comparable to the value related to 
the ionised reflection log Reflection = 3.0 + 0.2, suggesting a link 
between them. 

If the ionised reflection is associated to the warm absorber 
(e.g. the inner walls of a conical outflow), the lack of spectral 
features associated with the accretion disc is a natural conse- 
quence thereof. If, instead, the ionised reflection occurs at the 
accretion disc, it cannot extend up to the marginally stable or- 
bit. The lack of the significant relativistic blurring of this model 
component requires the disc to be truncated (inner disc radius 
Rin > 60 R g ). This idea is also supported by the low ratio of 
the reflection radiation to the primary one, R < 0.2, and also by 
the relatively low mass-accretion rate M « O.lMgdd determined 
from the source luminosity. 
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